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SECTION  I 
INTRODUCTION 

A  method  often  used  for  the  solution  of  initial-boundary  value  problems  in  structural 
mechanics  has  been  the  Finite  Element  Method  (References  l.and  2).  Many  applications  of  the 
Finite  Element  Method  in  linear  elastostatics  have  been  derived  from  the  Theorem  of  Minimum 
Potential  Energy  or  its  equivalent,  the  Principle  of  Virtual  Work  (Displacements),  Con¬ 
sequently,  the  set  of  primary  field  variables  is  the  displacement  field;  strains  and  stresses 
are  then  expressed  as  functions  of  the  displacements  through  the  field  equations. 

It  has  been  assumed  that  the  Finite  Element  Method  yields  displacements  and  stresses 
which  converge  as  the  mesh  size  is  decreased  (References  3,  4,  and  5).  Smooth  convergence 
of  the  displacement  field  has  been  widely  reported  for  many  applications  (Reference  2  a),  some 
violating  certain  of  the  compatibility  conditions  (Reference  1  a).  On  the  other  hand,  smooth 
convergence  of  the  stress  field  has  not  always  been  observed;  poor  convergence  and  spatial 
oscillations  are  found  (Reference  6).  Several  authors  have  given  stress  “averaging”  techniques 
(References  7  and  8).  For  many  engineering  purposes  it  must  be  recognized  that  a  knowledge 
of  the  stress  distribution  even  in  localized  regions,  is  the  desired  result  and  the  displace¬ 
ments  are  of  secondary  interest. 

Some  of  these  difficulties  can  be  circumvented  by  mesh  resolution  and  adoption  of  averaging 
techniques  mentioned  above  in  localized  regions  with  high  gradients.  However,  for  some  prob¬ 
lems  this  is  impossible  or  impracticable  and  other  alternatives  must  be  sought.  One  alter¬ 
native  is  to  employ  a  higher  order  displacement  expansion.  Investigators  have  pursued  this 
approach  with  encouraging  results  (Reference  9  and  10);  however,  the  stress  vector  remains 
discontinuous  over  contiguous  element  interfaces,  giving  a  stress  distribution  with  a  histograph 
appearance. 

Another  approach  is  to  obtain  the  stress  directly  as  a  primary  variable  through  an  applica¬ 
tion  of  the  Theorem  of  Minimum  Complementary  Energy.  Investigators  have  had  limited 
success  with  the  method,  one  difficulty  being  to  establish  a  stress  distribution  identically 
satisfying  the  equilibrium  equations  and  expressible  in  terms  of  a  convenient  set  of  global 
generalized  coordinates.  At  present  the  method  has  not  been  applied  with  the  same  degree  of 
generality  as  the  displacement  method  (References  2b,  11,  12  and  13).  It  is  surprising  that 
applications  of  mixed  variational  theorems  have  not  been  explored  in  this  regard  since  they 
offer  the  generality  of  the  displacement  method  and  retain  the  stress  field  as  a  primary 
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variable.  To  this  end  thepresentpaperintroduces  a  Finite  Element  application  of  the  Hellinger- 
Reissner  Variational  Theorem, 

It  is  shown  that  expansions  for  the  displacement  and  stress  fields  can  be  made  independently 
in  Finite  Element  analyses  and  that  continuity  of  the  stress  vector  can  be  maintained  across 
contiguous  element  interfaces.  Stresses  from  a  mixed  variational  theorem  are  shown  to  be 
improved  over  those  of  comparable  displacement  models.  The  limitation  principle  (Reference 
2b)  is  discussed  in  this  connection  and  is  shown  to  limit  the  scope  of  mixed  variational 
theorems. 

A  stiffness  matrix  for  plane  and  axisymmetric  solids  is  discussed  and  several  examples 
are  presented  with  comparisons  made  to  a  frequently  used  displacement  model  (Reference  17). 
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SECTION  n 

HELLINGER-REISSNER  THEOREM 


The  functional  in  the  theorem  is  expressed  in  terms  of  the  displacement  and  stress  field* 
(Reference  14) 
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The  theorem  asserts  that  of  all  admissible**  displacement  and  stress  states,  those  states  that 
satisfy  the  stress  equations  of  equilibrium,  the  strain  stress  equations  and  the  boundary 
conditions  render  the  functional  stationary 

Sv(u,  r  )  =  0  (3) 


Unlike  the  minimum  principles  the  stationary  value  cannot  be  shown  to  be  an  extremuum  value 
of  the  functional  (Reference  14). 


The  form  of  the  functional  is  such  that  separate  expansions  are  possible  without  the 
restrictions  necessary  in  the  minimum  principles.  The  stress  field  need  not  satisfy  the  equil¬ 
ibrium  or  stress-strain-displacement  equations  in  R,  nor  need  either  the  displacement  or 
stress  fields  satisfy  the  prescribed  boundary  conditions.  The  Euler  Equations  are  the  stress 
equations  of  equilibrium,  the  strain-stress  law  and  the  boundary  conditions.  Since  the  stationary 
value  of  the  functional  contains  the  constitutive  equation  as  an  Euler  Equation,  the  stress 
field  obtained  through  the  stress- strain-displacement  relations  may  not  agree  with  the  stress 
field  given  by  the  stationary  value. 


Standard  tensor  notation  is  used  and  a  cartesian  reference  frame  is  assumed;  comma  denotes 
differentiation,  repeated  indices  are  summed,  and  parentheses  denote  the  symmetric  part 
of  a  tensor.  Interpretation  of  the  symbols  is  given  in  the  Appendix. 


**  Admissibility  and  convergence  questions  are  beyond  the  scope  of  the  present  work  (Reference 
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If  independent  expansions  are  used  for  both  fields,  where  components  of  the  displacement 
and  stress  fields  are  adjoined  across  contiguous  element  boundaries  using  a  system  of  common 
displacement  and  stress  nodal  points,  then  a  direct  solution  for  the  primary  variables  should 
differ  from  corresponding  applications  of  the  minimum  principles.  However,  the  limitation 
principle  requires  that  the  stress  variables  be  left  free  {Reference  2b),  The  notion  of  con¬ 
nectivity  in  the  stress  variables  is  central  to  this  concept.  To  illustrate  this,  consider  the 
following*: 

ui  I(£j  (*p)o" 

r..  (xp)bn  (5) 


where  the  a’ s  and  b’s  are  the  generalized  coordinates  for  displacement  and  stress,  respectively. 
Suppose  there  exists  a  set  of  scalars  /3  ™  such  that 
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Then  if  there  is  no  interelement  connectivity  in  the  stress  variables,  Euler  Equation  of  the 
functional  Is  identically  satisfied.  The  stationary  value  yields,  in  part, 
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using  Equation  4  in  Equation  1,  and  assuming 
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Equation  7  becomes 
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*  The  superscripts  refer  to  generalized  coordinates  and  are  treated  as  ordinary  indices. 
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Equations  6  and  10  imply  that  the  stress  is  computed  from  the  stress-strain-displacement 
relations  and 


thus 


r.  u  . . 

i)  («,j) 
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in  which  case  Equation  12  and  Equation  8  imply  the  functional  in  Equation  1  can  be  reduced  to 

V(u)=f-L-u..c'  u  dv-f  P.  u.  da.  (13) 

2  ,,j  ,,kl  k,l  J  it 

r 

Equation  13  is  the  functional  in  the  Theorem  of  Minimum  Potential  Energy. 


The  above  is  a  demonstration  of  the  limitation  principle  which  is  seen  to  hold  in  the 
restrictive  case  where  the  stress  variables  are  not  connected  from  element  to  element. 
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SECTION  m 

A  FINITE  ELEMENT  MODEL 


Following  the  approach  in  (Reference  16)  a  linear  expansion  is  taken  for  both  the  stress 
and  displacement,  and  a  system  of  generalized  coordinates  is  established  at  the  vertices  of 
a  triangular  element  for  a  plane  stress  or  axisymmetric  solid.  Matrix  notation  is  employed 
for  convenience.  Displacement  boundary  conditions  are  to  be  satisfied  in  the  usual  fashion. 
Derivation  of  an  element  stiffness  matrix  for  a  two-dimensional  plane  stress  solid  is  given 
for  a  linear  triangular  element.  A  more  complete  exposition  for  both  plane  and  axisymmetric 
solids  is  given  in  (Reference  16). 


For  the  stated  application,  the  functional  in  Equation  1  can  be  expressed 

V(u  ,r)  =  J  (<r>|Vu|  -W(r)  -<  u>  jf  j )  dv  -  J  <u>|p"|  do  (14) 
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and  the  brackets  <  >  and  j  j  are  used  to  represent  row  and  column  matrices,  respectively. 

The  coordinate  functions  expansions  are 
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and  the  displacement  gradients  are 
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The  generalized  coordinates  are  expressed  in  terms  of  the  values  of  the  physical  components 
of  the  displacement  and  stress  fields  at  the  element  vertices,  (Figure  1) 


where 


Thus 


and 
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The  functional  becomes 


where 
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The  stationary  value  yields  the  governing  equations 
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The  linear  triangular  element  is  used  to  form  the  quadrilateral  element  shown  in  Figure  1. 
Condensation  is  then  performed  on  the  stiffness  matrix  in  order  to  reduce  the  bandwidth  of 
the  governing  Equations  9  and  17.  However,  in  the  case  of  the  present  mixed  model  this  cannot 
be  effected  directly  due  to  a  singularity  found  during  the  Gaussian  elimination.  What  is  done 


T*xl*  Tyyi’rxy^ 

LINEAR  TRIANGLE  FOUR-TRIANGLE  QUADRILATERAL 


Figure  1,  Four-Triangle  Quadrilateral 
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is  to  constrain  the  displacements  at  the  center  node  to  be  the  average  of  the  four  exterior 
nodes 


us  =  a  (ui  u2  +  u3  +  %  ) 


v5  =  *4  (vt  +  v2  +  v3  4  V4  > 


where 
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(32) 
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(34) 


The  constrained  stiffness  matrix  becomes 


[Srj]  8  [Sru][*]  (35) 

v 

Standard  techniques  are  used  to  solve  the  resulting  system  of  linear  algebraic  Equation 
30.  Although  Equation  30  is  not  a  positive  definite  system,  no  difficulty  is  encountered  in 
Gaussian  elimination  if  nodal  point  stress  are  eliminated  before  the  displacements. 
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SECTION  IV 
EXAMPLES 


To  illustrate  the  method,  two  example  problems  are  given  with  direct  solutions  obtained 
using  the  mixed  model  and  a  comparable  displacement  model  Reference  17, 


A  boundary  value  problem  with  an  applied  linear  stress  distribution  is  shown  in  Figure  2. 
The  problem  is  equivalent  to  combined  extension  and  flexure  in  elementary  beam  theory.  The 
two-dimensional  plane  stress  solution  is 


y  v'i 

€  =  -=■  ,  €  = - —  ,  €  =0 

xx  E  yy  E  xy 

,  s  /  *  1  .  2  ,  2 

u '  x  ,  y )  =  — —  ,  v  (x  ,  y  )=-——•(  X  -H  j/  y  ) 

“  ^  U 


(36) 


The  linear  stress  distribution  is  contained  within  the  expansion  for  the  stress  field  in  Equation 
16,  but  the  displacement  field  is  parabolic. 


Figure  2  shows  a  plot  of  the  stress  field  along  the  top  fiber  of  the  beam  and  through  the 
depth  at  the  edge  of  applied  traction.  The  mixed  model  gives  accurate  values  for  the  stress 
field  including  those  points  on  the  boundary.  The  displacement  model  gives  values  for  the 
stress  only  at  the  element  centroid.  The  mixed  model  also  predicts  the  deformed  profile 
with  considerable  accuracy. 


Figure  3  illustrates  a  finite  length  hollow  cylinder,  fixed  at  the  outer  radius  and  acted  upon 
by  an  internal  pressure.  The  geometry  and  physical  properties  are 


a  =  4,  b  =  6,  i  --  5.5 
Y  =  I  ,  v  --  0  A 


(37) 


A  plot  of  the  normal  stress  distribution  along  the  fixed  boundary  is  given  in  Figure  3  for  the 
two  Finite  Element  models.  The  stress  singularity  is  seen  to  be  more  severe  in  the  mixed 
model,  the  displacement  model  suppressingmuch  of  the  gradient  in  the  vicinity  of  the  singular¬ 
ity. 
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[17] 


Figure  2.  A  Linear  Stress  Problem 
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Trr(D)  /  p 


Figure  3.  Pressurized  Finite  Cylinder 
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SECTION  V 
CONCLUSIONS 


For  the  examples  considered  the  present  application  of  the  Hellinger-Reissner  Theorem 
gives  a  more  accurate  field  description  of  both  displacement  and  stress  than  do  existing 
applications  of  displacement  models.  The  mixed  model  yields  solutions  in  which  all 
components  of  the  stress  tensor  are  continuous  from  element  to  element,  eliminating  the 
histograph  distribution  often  found  in  displacement  models.  The  displacement  vector  ob¬ 
tained  from  the  mixed  model  is  considerably  more  accurate  than  that  given  by  comparable 
displacement  models.  Results  for  the  stress  field  are  superior,  but  tend  to  be  dependent  on 
material  properties  and  mesh  configuration.  In  part,  this  can  be  attributed  to  identically 
satisfying  the  displacement  boundary  conditions.  Since  the  stress  distribution  is  the  desired 
result,  it  would  be  advantageous  to  bias  the  model  in  this  direction.  An  alternate  form  of  the 
functional  which  accomplishes  this  is  currently  being  studied.* 

The  mixed  model  is  shown  to  be  highly  effective  in  capturing  steep  stress  or  displacement 
gradients  that  can  occur  near  singularities  in  boundary  value  problems.  These  are  often  not 
predicted  by  displacement  models.  The  mixed  model  is  most  effectively  used  for  this  class  of 
problems  since  an  excess  of  generalized  coordinates  can  require  considerable  computational 
effort  for  a  general  problem. 

The  present  effort  should  be  regarded  as  introductory,  and  considerable  additional 
research  is  indicated.  Alternate  forms  should  be  considered.  The  number  of  generalized 
coordinates  should  be  reduced  by  employing  other  coordinate  functions  and  element  con¬ 
figurations.  Herrmann  has  shown  that  this  can  be  done  for  plate  bending  (Reference  lb),  and 
Prager  (Reference  18),  has  extended  and  systematized  the  relevant  principles.  What  has  been 
clearly  demonstrated  here  is  that  Finite  Element  Methods  can  be  successfully  used  with  mixed 
variational  theorems  for  direct  solutions  to  boundary  value  problems  in  mechanics. 


*This  form  has  been  suggested  by  S,  Pawsey  in  private  communications. 
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